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ABSTRACT 


A major goal of the hydraulic ram survivability program 
for aircraft fuel tanks is the development of analytical- 
numerical tools for the accurate prediction of damage to the 
tank due to a ballistic projectile. This report presents a 
method for predicting the amount of cracking of a penetrated 
tank wall due to the penetrating projectile and the hydraulic 
ram loading. The method uses computer codes to predict the 
fluid pressure on the wall and the stresses in the wall. The 
stresses are compared with empirical data on the fracture of 
thin cracked plates to obtain a prediction of the final crack 
length. A comparison of predicted cracks with actual cracks 
that occurred in plates tested by Naval Weapons Center, China 
Lake, California is included. Good correlation is obtained 
when the magnitude of the predicted strains is adjusted to 
agree with the magnitude of the measured strains by a correc- 


teen factor. 
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I. INTRODUCTION 


A. FUEL TANK DAMAGE 

Studies of aircraft battle damage reports from the Viet 
Nam War have pointed out the vulnerability of sophisticated 
aircraft to small arms ground fire and to surface-to-air 
Missiles. Since aircraft fuel cells have the largest sur- 
face area and volume of all the vulnerable wirenart components, 
and in light of increased aircraft costs and continued decreases 
in the national defense budget, aircraft fuel cell surviv- 
ability analysis and design are absolutely necessary. 

Since projectiles represent one of the greatest threats 
to aircraft fuel tanks, projectile-fuel tank interaction is 
of particular interest. A projectile that penetrates a 
fluid-containing cell causes damage many times greater than 
a projectile penetrating an empty cell. This interaction of 
a projectile with a fluid-containing cell is called the 
hydraulic ram phenomenon. 

The hydraulic ram phenomenon can be divided into five 
phases which are named after the dominant mechanism of pro- 
jectile energy dissipation. The phases are: 

1. “Penetration phase" ... projectile energy is lost 

in "punching" a hole in the entry wall. 

2. "Shock phase" ... projectile protrusion into the 

fluid drives a hemispherical shock wave into the 


tank liquid. 
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3. "Drag phase" ... projectile energy is dissipated 
by form drag as it passes through the fluid, and a 
cylindrical air- and vapor-filled cavity is formed. 

4, “Exit phase" ... projectile punches a hole in the 

exit wall if first three phases have not dissipated 
all of the projectile energy. 

>. Cavity oscillation phase"). results trometne 

growth and eventual collapse of the drag phase 
cavity. 
The interrelationships of these phases are shown in Fig. 1, 
and more detailed explanations of each phase may be found in 
feme. 1, 2, 3 and 4.§ It should be noted that there is no 
precise beginning or end of any particular phase. 

The complete penetration of a fluid-containing cell by a 
ballistic projectile can have a catastrophic effect on the cell 
exit wall, as can be seen in Figs. 2 and 3. However, when 
the projectile velocity is lowered slightly comparatively small 
damage may occur as can be seen in Figs. 4 and 5. Although 
the shock phase and the cavity oscillation phase can affect 
the exit wall response, the drag phase and the exit phase are 
usually the most important when considering damage to exit 
walls. During these two phases, the exit wall is loaded by 
pressures due to both the penetrating projectile itself and 
the fluid motion toward the wall caused by the earlier shock 
and drag phases. These pressures cause the relatively large 


and small damage shown in Figs. 2-5. The ability to predict 


13 





y 


CESS erg 
HUTMOTTOI osTNg sanssoidg 


\ 
\ 
N 


ARTARPD FO osde{[Too 


YURL PeT{Ita ptnbtT 


OTTROSLOA 


Hutjtquny Aq poyerous 
SABM OANSSerg pue ARTA 


[ @_/ NONSWONHHd WVU OIINWECAH 





T dandTd 





yUeL PSTTTA prnbt7T 


RoVedwu~ sTTRzOeLlOIg 

Aq pesneg oansserg ySTH 
5ButAueduosoy pue oAeM Yyoous 
yueL PeTTtda prnbtyT 





14 











FIGURE 2. NWC 6-HR-3 Series Test Plate 
Thickness 0.063 inch 
Projectile velocity 2358 pS 





FIGURE 3. NWC 6-HR-14 Series Test Plate 
Thickness 0.063 inch 
Projectile velocity 2212 fps 
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FIGURE 4. NWC 6-HR-8 Series Test Plate 
Thickness 0.063 inch 
Projectile velocity 2118 fps 








FIGURE 5. NWC 6-HR-15 Series Test Plate 
Thickness 0.063 inch 
Projectile velocity Bie, pS 
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the amount of damage to the exit wall is important in fuel 
cell survivability analysis and design because: 

1. a prediction of the time required for the fuel in 
a tank to empty can be made. 

2. designers can predict the amount of fuel dumped into 
an engine inlet from an integral fuselage-intake fuel 
tank. 

3. designers can estimate the amount of fuel atomization 


TtAGOcelinogwanduatne associated fire hazard. 


B. PREVIOUS HYDRAULIC RAM AND STRUCTURAL RESPONSE STUDIES 


Due to the increased emphasis that has recently been placed 
on aircraft survivability, a large amount of research effort 
has been devoted to the development of an understanding of the 
hydraulic ram effect and the associated fluid-structure inter- 
action. The Hydraulic Ram Project conducted at the Naval 
feapoens Center (NWC), China Lake, California, is of particular 
interest here. Phase One of that project resulted in a digital 
computer code, Ref. 5, developed by Lundstrom and Fung that 
predicts the drag phase fluid pressure throughout a rectangu- 
lar tank due to ballistic penetrators. The code was based 
upon an analytical model, Ref. 2, and the empirical data of 
Ref. 6, which was also obtained during Phase One. 

Phase Two of the NWC Hydraulic Ram Project was concerned 
with the description of the fluid-structure interaction during 
the loading of the fuel tank exit wall due to the hydraulic 


ram pressure. In support of this goal, detailed measurements 


Jey 





of the strains at several locations on the tank exit wall 
were obtained from an extensive series of v7 Mime ek balls = 
tic penetration tests on rectangular, fluid-filled tanks. A 
detailed description of the test Set-up and the results are 
given in Ref. 7. 

An analytical and experimental hydraulic ram program has 
been conducted at the Naval Postgraduate School (NPS) in 
conjunction with the NWC program for approximately three 
years. As part of this program a ballistic range, consisting 
of a 22 caliber rifle and a fluid-containing cubic tank, was 
built and used to obtain fluid pressure and entry wall strain 
measurements for various projectile sizes and energy levels. 
The results of the various experiments conducted are given in 
Refs. 8-18. 

The goal of the analytical phase of the NPS program has 
been to develop methods of analysis that will accurately pre- 
dict the structural response of the tank walls to loadings 
caused by hydraulic ram. The results of this effort are given 
in Refs. ll and 18-22. Two structural analysis digital com- 
puter codes, BR-1 and SATANS, have been modified by Ball to 
predict the stresses and strains in the walls of a fluid- 
containing tank as a function of time. It was necessary to 
modify the two computer codes in order to account for the 
fluid-structure interaction and cavitation. One-dimensional 
piston theory was used to approximate the fluid-structure 
interaction. Information on the modified codes, BR-1HR for 


rectangular walls and SATANS for circular walls, and their 


eS) 





use 1s given in Ref. 23 and 24. Both BR-1HR and SATANS use 
the NWC Hydraulic Ram Pressure code described in Ref. 5 to 
determine the pressure applied to the tank wall. A compari- 
son of the exit wall strains predicted by SATANS with the 


strains measured by NWC is given in Ref. 25. 


C. EXIT WALL FRACTURE STUDY 


This report presents the results of research directed at 
understanding and predicting the fracture of fuel tank exit 
walls caused by the hydraulic ram phenomenon. The information 
in Section II begins with a description of several common types 
of fractures. Section II continues with a few of the current 
theories of fracture that were considered as possible theories 
for predicting the extent of fracture in the exit walls. In 
section III is given a brief description of the hydraulic ram 
exit wall tests conducted at NWC for three different velocity 
levels. In Section IV, the fractured plates, strain gauge 
readings, and high speed motion pictures from the NWC tests 
are examined in order to determine the fracture process and 
to guide the development of a fracture prediction method. In 
Section V, the proposed method of fracture prediction is 
Mecribed, and a comparison of the predicted fracture length 
to the actual fracture length that aeeurned in the NWC low 
velocity tests is presented. Section VI contains the conclu- 


sions of this study and recommended areas for further study. 
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II. FRACTURE AND THEORIES OF FRACTURE 


An understanding of fracture and the theories of fracture 
is important when predicting the damage to the fuel cell exit 
wall. When stress causes a material to break into two or 
more parts the material is said to have fractured. Although 
fracture is a relatively common occurrence, its true nature 
is still not completely understood. One of the most convenient 
ways to describe fracture and the theories of fracture centers 
See tie Size of the fracture region of interest. There are 
usually three size levels of concern; atomistic level (dimen- 


8 


Sions on the order of 10 ~ in.), microscopic level (dimensions 


4) 


on the order of 5x10 “~ in.), and macroscopic level (dimensions 
on the order of 107+ in. or greater). However, as pointed out 
by Drucker / Ref. 26 /, the atomistic and microscopic levels 
are better suited to describing fracture than predicting it. 
Also, because the accuracy of the predicted stress-strain 
histories at these two levels is debatable, all three levels 


will be used to describe fracture, but only the macroscopic 


theories of fracture will be discussed. 


pee LYPES OF FRACTURE 

At the atomistic level individual atoms and molecules are 
the center of attention. At this level, fracture consists of 
two types; cleavage and shear. Cleavage is defined as the 
breaking of atomic bonds normal to the fracture plane, and 
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shear is defined as the process of shearing the bonds across 
the fracture plane. 

At the microscopic and macroscopic eee the discussion 
becomes more complicated. The microscopic level centers on 
the individual crystals of material, and the macroscopic level 
centers on the total specimen. In both cases there are three 
classes of fracture; cleavage, shear, and grain boundary. The 
first two can be described as transcrystalline and the third 
can be described as intercrystalline. Grain boundary or inter- 
crystalline fracture usually occurs under a static loading 
condition that leads to creep fracture and failure due to 
stress corrosion. Because this loading condition and these 
types of failure are not found in the hydraulic ram process, 
intercrystalline fracture will not be considered further. 

Microstructurally, cleavage fracture can be distinguished 
from shear fracture by a consideration of the characteristic 
Lattice plane on which the fracture occurs. To illustrate 
which plane is associated with the two types of fracture, 
consider the unit cell of a body centered cubic crystal such 
as found in iron or steel, shown in Fig. 6. The cleavage 
plane / 100_/ of Fig. 6(a) locates the cleavage fracture while 
it is the slip plane / 110_/ of Fig. 6(b) on which the shear 
fracture operates. Table I gives the cleavage and shear 
planes for various crystalline structures and materials. As 


can be seen in Table I, aluminum (Al) does not have a cleavage 


plane. 
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(100) cleavage plane GbE) eilije joleiae 


PIGURE 6. Characteristic transcrystalline fracture planes 
in iron or steel 7297 


TABLE I 


Cleavage and Shear Planes for 
Various Crystal Structures and Materials / 28 7 





Crystal Cleavage Ease, 
Structure Example Plane Shear Planes 
Livia, K,V,Cr, 
BCC Fe,most steels, [100] [112], [110] 
Min Cb, Moy, ra 
Cu,Ag,Au,AlNi, | 
FCC brass,300 series None {111] 


stainless steels 
Be, Mg, 2n, 


HCP Sn ari., Oe [1000] [1122], (kone, 
Cai, graphite [1000} 


Ze 





Fractography is the macroscopic and/or microscopic exami- 
nation of the fracture surface and is used to provide infor- 
mation on the micromechanism of fracture. On a macroscopic 
level the general appearance of the cleavage fracture surface 
is bright, granular, and crystalline while the shear fracture 
is dull gray, fibrous, and silky. On the microscopic level 
the examination of fracture surfaces can be accomplished with 
elther a light microscope or scanning electron microscope (SEM). 
The interested reader is directed to Ref. 29 for a complete 
discussion of microfractography and its interpretation. 

Figure 7 shows the general appearance of cleavage and shear 
fractures on the macroscopic and microscopic levels under uni- 
axial tensile loading. Macroscopically, besides the general 
appearance already discussed, cleavage fractures appear flat 
or square as shown in Fig. ViajOm the other hands) shea 
fractures appear either slanted or chisel-pointed as shown in 
Bere, 7(b). 

Another common way of distinguishing different fracture 
types is based on the amount of plastic flow involved in the 
fracture process. Ductile fracture is characterized as having 
a considerable amount of plastic flow associated with it. 
Brittle fracture, in contrast, is characterized as having 
very little or no plastic flow. Because these two terms are 
not consistently associated with any single one of the terms 
already discussed, they greatly confuse the issue. Table Il 
may be helpful in elucidating the terminology, if it is 


remembered that the two right-hand terms of the third row 
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Macroscopic level : s 
Microscopic level 






oO 
| Continuous 
cleavag 
(a) 7 a 
Discontinuous 
cleavage 
Cleavage 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 





(D) 1 
Slant Chisel 
fracture: Pola: | shear rupture 
fracture 
Shear or shear rupture 
FIGURE 7. Comparison of Cleavage and Shear Fractures 


at Macroscopic Level and Microscopic Level 
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aba 


Transcrystalline fracture and its terminology / 27_/ 


Basis of classification Types of transcrystalline fracture 
Crystallography shear cleavage 
Fracture appearance fibrous granular 
Strain or strain energy at ductile brittle 

imeac Lure 


are often used either interchangeably with the term of the 
first row in the same column or as adjectives describing 
either term of the first row / 27_/. 

Since fuel cell walls are often thin sheets of material, 
the fracture of sheet material in tension is of interest. 
Tetelman and McEvily / Ref. 28_/ and Ryder and Smale / Ref. 
29_/ describe the two macroscopically different types of 
behavior illustrated in Fig. 8. Both types of behavior are 
characteristic of shear fracture. Figure 8(A) shows a 
relatively brittle process that would be expected in thick 
plates and round tensile bars. Figure 8(B), on the other 
hand, shows a relatively ductile process which would be 
expected in thin sheets, particularly those sheets made of 


the Duralumins. 


B. MACROSCOPIC FRACTURE THEORIES 

Many theories of fracture have been developed at the 
Macroscopic level. The various theories range from simple 
to complex in nature. Some assume an unflawed specimen, and 
others do not. Only a few of the more widely accepted the- 


orlies will be described. 
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& Wa Direc wore 

9 Crack Growth 
y - 
(A) 


(B) 


FIGURE 8. Two Modes of Tensile Fracture Observed in 
sheet Specimens: (A) Brittle Normal Fracture; 
(B) Ductile 45~ Fracture LEO 





The first three theories are similar in that they assume 
a specimen with no defects and use either the tensile (or 
compression) ultimate strength, De , OF tne) 1) tama temchic ar 
strength of the material as the criterion for fracture. The 
Simplest of these macroscopic theories is the maximum tensile 
stress, or Rankine, theory. This theory would predict fracture 
if the maximum principal stress of a given specimen was greater 
than or equal to the tensile (or compressive) ultimate strength. 
The theory is generally considered reasonable only when failure 
occurs by cleavage fracture. 

The maximum shear stress, or Tresca, theory proposes that 
fracture will occur when the maximum shear stress exceeds the 
ultimate shear strength of the material. When used to predict 
the fracture of polycrystalline materials, such as aluminun, 
this theory Baty also requires the definition of a criterion 
mductility. 

The last theory of this type, closely related to the maxi- 
mum shear stress theory, is the octahedral shear stress theory. 
This theory is also known as the maximum distortional energy 
theory or the von Mises theory. In terms of the principal 


stresses Orne and O, this theory may be expressed as 
Ze 2 Za 2 
aO> 07) “= (O> Ose Os ©. = 267: 


This theory is generally used for materials that exhibit 
large amounts of plasticity or flow before fracture. The 


fracture criteria for these three theories in plane stress 
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are compared in Fig. 9 by plotting them on a graph of biaxial 


stress OF and Oy. 


Maximum Stress 









Maximum 
shear 
(TRESCA) 


Octahedral 
shear 
(vonMises) 


FIGURE 9. Comparison of Strength Theories for biaxial 
; stress 
lf a defect, such as a crack or hole, exists in the speci- 


men, the fracture analysis must determine whether or not the 


crack will propagate under a given load. This problem is 


usually treated using the modifications of the Griffith 


crack theory proposed by Irwin / Ref. 31/7 and Orowan / Ref. 32_/. 
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These modifications led to a concept of fracture instability 
based on a stress-intensity factor kK aC sphie cOOcmow asa 
tionary or slowly growing crack or defect in a body. The 
elastic confinement of the material governs the degree of 
plain strain at the tip of the crack. Therefore the critical 
stress intensity, Ky is assumed to be a material property. 
The theory proposes that when the loading produces a stress 
intensity equal to the critical one, the oraee propagates 
rapidly over the cross section. The critical stress-intensity 
factor can be related to the strain-energy release rate since, 
at the instability of crack extension, the energy release rate 
is equal to the rate of change of strain energy with respect 
to crack area when the displacements of the externally applied 
loads are held constant during the incremental crack extension. 
Both the strain-energy release rate and the stress-intensity 
factor are functions of loading, structure geometry and size, 
and crack orientation. 

The stress-intensity factor also varies with the relative 
motion of the crack surfaces,or what is called the modes of 
fracture. There are three basic modes or combinations of 
them. Mode I, the opening mode, is characterized by the crack 
surfaces separating, perpendicular to the maximum tensile 
stress. Mode II, the shear mode, is characterized by the 
crack surfaces shearing in a direction perpendicular to the 
crack front. Mode III, the tearing mode, is characterized by 
the crack surfaces shearing in a direction parallel to the 


crack front. For example, the stress-intensity factor 
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Ky =O (Ta)? 


is for an infinite plate under a uniaxial average tension 
stressO with a crack length of 2a perpendicular to the 

direction of applied stress. The subscript I denotes the 
opening mode. Fast crack extension occurs when Ky = ee 


Thus, the critical stress for fracture Ox is 


Be 
O, 7 Ky,( a) ° 


Note that Ox varies inversely with the square root of the 
crack length. Therefore, if the crack is sufficiently long, 
fracture will occur when the average stress in the plate is 
below the yield strength. 

Modifications to include plasticity around the crack 
tip make it possible to handle materials that exhibit semi- 
brittle or ductile fractures. The reader is directed to 
Ref. 33 for a more complete discussion of this theory and its 
various modifications. 

The fracture criteria discussed in this Section were 
each considered for use in the crack prediction method de- 
scribed in Section V. Those theories that pertain to unflawed 
specimens such as the maximum shear stress theory or the 
octahedral shear stress theory were disregarded because they 
are generally only considered relevant for unflawed specimens. 
The modified Griffith crack theory described was not used 


because the selection of the appropriate modifications to 
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mame theory for the test plate feometry Considerations, and 
the selection of an appropriate plane stress critical stress 
intensity, Kos to use, was not considered feasible at the 
present time. Since the theoretical fracture criteria did 
not appear useable for the reasons mentioned, a method based 
on empirical data was formulated. The method is presented in 


section V. 
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IJ]I. THE NWC EXIT WALL TEST PROGRAM 





A. EXPERIMENTAL SET-UP 


During the last three years the Naval Weapons Center (NWC), 
China Lake, California, has conducted experimental work on the 
hydraulic ram phenomenon which is directly related to this 
study. In this section the relevant part of this work will 
be briefly discussed. 

A series of tests was conducted at NWC to measure the 
strains at four locations on the exit wall of a simulated 
aircraft fuel tank subjected to hydraulic ram effects. The 
tests were conducted at NWC, using 12.7 mm API ballistic 
projectiles fired at a fluid filled, rectangular tank. The 
tank is 60 inches wide by 60 inches high, with 22 inches 
between the entry and exit walls, and has an open top. A 
schema of the tank is shown in Fig. 10. The central portion 
of the entry wall is a stretched rubber membrane and the central 
portion of the exit wall is a 20-inch by 20-inch by 0.125-inch 
test plate of 2024-T3 aluminum clamped to a 0.25-inch steel 
wall. The exit wall with a damaged test plate is shown in 
Fig. 11. The edge clamping of the test plates was obtained 
by compression (friction) between two rubber gaskets around 
the outer perimeter of the plate. A pin through the plate at 
two diametrically opposed corners provided additional edge 


support. A schema of the test set-up is shown in Fig. 12. 
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FIGURE 10. Schema of Test Tank 
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FIGURE 11. Test Plate Mounted in Exit Wall 
of Test Tank 
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FIGURE 12. Schema of Test Set Up 
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During this series of tests fourteen 12.7mm API projectiles 
were fired at the tank with velocities between 1300 fps and 
2800 fps. The strains on the wet and dry surfaces of the 
exit wall were measured at the four locations shown in Fig. 13 
for each shot. The strains at each location were measured in 
the two directions indicated by the short lines in Fig. 13. 
Pressures were measured at five locations along the projectile 
trajectory inside the tank. High-speed motion pictures of the 
exit wall response were taken looking at the dry side of the 


exit wall. 


B. EXIT WALL DAMAGE 


Table III lists the projectile velocity and the extent of 
damage for each shot, and Figs. 14-26 show the damage looking 
at the dry side of each test plate. In general, the low veloc- 
ity projectiles caused limited bulging of the exit wall with 
short cracks approximately three inches in length emanating 
from the penetration point. The high velocity projectiles, 
on the other hand, caused severe tearing of the wall with 
several cracks emanating from the exit hole and running to 
the edges of the plate. A more detailed description of the 
experimental set-up and a presentation of all exit wall strain 
data, is given in Ref. 7. Selected ieee from the NWC 
tests will be presented in Section V where a comparison be- 
tween the experimental results and analytical results is 


made. 
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Velocity (fps) 
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Projectile Velocity and Extent of Exit Wall Damage 


Damage 
Cigar-Shaped hole 3 in. long 
Sia MinGhoetteamanere 1.5 ein. 1n 
dia. with one 5 in. crack and 
one 1 in. crack 
One 1l in. and one 9 in. cracks 
Cigar-shaped hole 3 in. long 
Cigar-shaped hole 3 in. long 
Circular hole 3 ineetnecaa. in- 
cluding petaling with one 10 in. 


Crack and one 8 in. -crack 


5 in. hole with several 5 in. 
cracks and petaling 


5 in. hole and several 5 in. 
cracks and petaling 


Bulged plate with 5 in. cigar- 
shaped hole and two 1.5 in. 
cracks 


One 3 in. crack and two 1 in. 
cracks 


Cigar-shaped hole 2 in. long and 
two 2 in. cracks 


Cigar-shaped hole 2.5 in. long 
and four i.5 1n. Cracks 


Not available 


Very small hole with one 3 in., 
One 2 In<,eanad One Iso in. Glracims 
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IV. FRACTOGRAPHIC EXAMINATION OF EXIT WALL 
FRACTURE SURFACE 

One of the major objectives of this study was to determine 
the type of fracture caused by the hydraulic ram pressure 
loading, and if possible to determine the significance of 
eu@ain rate effects on the fracture process. Consequently, 
a fractographic examination of the fracture surfaces of the 
NWC test plates was carried out. The photographs of the 
damaged plates, Figs. 14-26, and the actual NWC 7-HR-7 test 
plate were examined macroscopically, and various pieces of the 
fracture surfaces from the 7-HR-7 test plate were examined 


with the NPS scanning electron microscope (SEM). 


A. MACROSCOPIC EXAMINATION 

Examination of the fractures shown in Figs. 14-26 and 
the 7-HR-7 test plate showed that the fracture surfaces are 
generally slanted at approximately 5? to the plate surfaces. 
As described in Section II and shown in Fig. 10(a), this 45° 
Slant of the fracture surface is characteristic of a relatively 
ductile fracture process of a thin plate in a plane stress 
condition. The Te slant of the fracture surface is also an 
indication of a shear fracture, as described in Section II. 
The identification of the fractures as shear fractures was 
further strengthened by the fiborous, dull gray appearance 


ee the 7-HR-7 test plate fracture suriaces. 


pe 





B. SCANNING ELECTRON MICROSCOPE EXAMINATION 

The examination of the funcwtee surfaces of the 7-HR-7 
test plate with the SEM was considerably more complicated 
than the macroscopic examination. Figure 27 shows the loca- 
tions on the 7-HR-7 test plate where the SEM specimens were 
taken. Figures 28-49 are the fractographs obtained using the 
SEM. As seen in Figs. 28-49, what appeared to be of uniform 
Character to the unaided eye becomes extremely varied in 
nature when magnified 500 to 1200 times. These microfractographs 
have been compared to the microfractographs of 2024-13, 2024-TH, 
and 7075-T6 aluminum shown in Ref. 29. The fractographs pre- 
sented in Ref. 29 were obtained from samples fractured under 
very different circumstances than those of the 7-HR series, 
and are of different materials in some cases. Nevertheless, 
they are very similar to some of the fractographs presented 
here. 

Figures 28-31 are SEM fractographs of specimen number 1. 
This specimen was obtained from location 1 shown in Fig. 27. 
This surface is where the projectile penetrated the plate. 
There were no fractographs in Ref. 29 similar to these 
four fractographs. Macroscopically this specimen appeared 
different from the others. The difference is probably due 
to the fact that the projectile penetrated the plate at this 
location. The surfaces shown in Figs. 28-31 appear to have 
been produced by a combination of fracture and scraping as the 


projectile penetrated the plate. 
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mr GURE Ze on of SEM Specimens on 7-HR-7 Test 
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FIGURE 28. SEM Fraetograph of Specimen #1 ...x620 





FIGURE 29. SEM Fractograph of Specimen #1 x600 
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FIGURE 31. 


SEM Fractograph of 
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FIGURE 32. SEM Fractograph of Specimen #2 x560 





FIGURE 33. SEM Fractograph of Specimen #2 x110 
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FIGURE 34. SEM Fractograph of Specimen #2 x 550 








FIGURE . SEM Fractograph of Specimen #2 x1100 
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FIGURE 36. SEM Fractograph of Specimen #2 x550 





FIGURE 37. SEM Fractograph of Specimen #2 x1100 
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FIGURE 40. SEM Fractograph of Specimen #4  x600 








FIGURE 41, SEM Fractograph of Specimen #4 x1250 
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FIGURE 42. SEM Fractograph of Specimen #4 X625 








FIGURE 43. SEM Fractograph of Specimen #5 x1200 
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FIGURE 44, SEM Fractograph of Specimen #5 x600 








FIGURE 45, SEM Fractograph of Specimen #5 x1100 








FIGURE 46. SEM Fractograph of Specimen #5 x550 








FIGURE 47. SEM Fractograph of Specimen #3 x1200 
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FIGURE 48. 


SEM Fractograph of Specimen #3 





FIGURE 49. 


SEM Fractograph of Specimen #3 
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Figures 32-46 are SEM fractographs of specimens 2 ail 
5, as indicated in the captions. These figures are similar 
enough to be considered together. These figures were compared 
with Figs. 4022, 4023, and 4028 of Ref. 29 and that portion 
of Ref. 29 entitled "Interpretation of Scanning-Electron- 
Microscope Fractographs." Most of the surfaces of the speci- 
mens are composed of dimples formed by microvoid coalescence 
with tear ridges. The dimples are indicative of a shear type 
fracture with the dimples and tear ridges also indicating a 
ductile process. Some of the large smooth areas, such as 
those in Figs. 35, 37, and 43, could possibly be cleavage or 
intercrystalline fractures, but these areas are believed to 
be the result of microvoid surfaces that have been stretched. 
The elimination of the possibility of cleavage fracture seems 
justified due to the lack of river patterns in the large smooth 
areas which are characteristic of cleavage fractures. The 
elimination of the possibility of intercrystalline fracture 
seems justified due to the lack of intergranular secondary 
fissures such as shown in Fig. 4028 of Ref. 29. 

Figures 47-49 are SEM fractographs of specimen 3 which, 
like specimen 1, came from very close to the projectile pene- 
tration point. The surface represented by these fractographs 
could not be characterized into a definite pattern, although 


there is an indication of some dimpling. 
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Sel LGH SPEED MOTTON PIChuRES 

High speed motion pictures taken of the 7-HR series tests 
showed lapsed time from projectile penetration of the plate 
until the cracks stopped increasing of between 0.625 and 2.0 
msec. The longer times were associated with the 7-HR high 
velocity, high damage shots 6-9, and the shorter times were 
associated with the 7-HR low velocity shots 10-14. The high 


speed films are discussed.in more detail in Section V. 
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V. EXIT WALL CRACK LENGTH ANALYSIS AND RESULTS 


The second major objective of this study was to develop a 
macroscopic method for predicting the extent of cracking in 
the exit wall of a fuel tank due to hydraulic ram. Shots OR 
11, 12, and 14 of the NWC 7-HR series tests were used in the 
development of the proposed prediction method because of their 
low velocity and the relatively minor damage they caused. The 
method of crack length prediction follows. 

Computer codes for the fluid pressure and the structural 
response are used to predict the stresses in the exit wall. 

A correction factor, determined as explained later in this 
section, is applied to the predicted stresses. The corrected 
stresses are plotted as a function of the distance from the 
projectile penetration point. Reference is made to empirical 
data to find the stress required to increase the crack length 
for various initial crack lengths. The stress required to 
increase the crack length is plotted as a function of the 
crack length on the same graph as the corrected stresses. 
Thus, the predicted crack length is that crack length at 
which the corrected stress first falls below the stress 
required for crack extension. 


A. ANALYSIS FOR MAXIMUM PREDICTED WALL STRAINS AND A 
COMPARISON WITH EXPERIMENTAL DATA 


The NWC digital computer program for the hydraulic ram 


wall pressures, as described in Ref. 5, and SATANS, as modified 
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in Ref. 24, have been used to predict the axi-symmetric 
strains in the exit wall on a radial line at one-half inch 
intervals. In the analysis, the projectile penetration point 
was assumed to be at the center of the plate, as shown in Fig. 
50. SATANS was used for the analysis of the rectangular wall 
instead of BR-1HR because its computer execution time is consid- 
erably less than that of BR-l1HR, which takes several hours of 
IBM 360/67 execution time for a typical solution using a coarse 
model. The parameters used in the two computer codes were 
those of the NWC 7-HR-14 test, which is similar enough to the 
7-HR-10, 11, and 12 shots to make the results comparable for 
any of the four tests. Piston theory was used to approximate 
the fluid-structure interaction. Cavitation at the exit wall 
was assumed if the net interface pressure between the fluid and 
the wall became negative (tension). The results of the com- 
puter analysis are compared with the experimentally measured 
Stearns in Ref. 25. 

According to Ref. 25, the magnitudes of the strains pre- 
dicted by the computer analysis do not agree very closely with 
the actual strains measured in the NWC 7-HR series experiments. 
Figures 51-56 are graphs of the actual strains measured by the 
strain gages located as shown in Fig. 13 and of the predicted 
strains at the point marked A in Fig. 13 and Fig. 50. Fig- 
ures 51, 52, and 55 are the strains, actual and predicted as 
indicated, in the y-direction. Figures 53, 54, and 56 are the 


strains, actual and predicted as indicated, in the x-direction. 
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FIGURE 53. Strains Measured by Gage 11 in 
7-HR-14 Shot 
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Although the magnitudes of the actual strains are much larger 
than the predicted strains, the basic shape of the correspond- 
ing strain histories are roughly similar. 

The maximum predicted strains at point A have been compared 
with the average of the maximum strains recorded on the corres- 
ponding strain gages at the two locations. The comparison 
showed that the x-direction predicted strains were low by 
a factor of 4.6 and the y-direction predicted strains were low 
by a factor of 2.4. One reason for the difference in correction 
meewors could be that SATANS analyzes a circular plate with 
the penetration point at the center, while the actual plate 
1s square with the penetration point located off center as 
shown in Fig. 50. In addition, the actual locations of the 
strain gages do not correspond exactly to the location of 
the predicted strains. Other possible reasons are that the 
method used to treat the fluid-structure interaction, the 
assumption of cavitation, and the hydraulic ram pressure 
analysis are not representative of the actual situation. A 
further discussion of the comparison is given in Ref. 25. 


The predicted strains were used to determine the stresses 


OF 


. and O, on the outside surface of the tank exit wall (dry), 


with © aligned circumferentially and Oy aligned radially. 
The maximum positive O, and OQ), stress predicted by SATANS at 
stations 2 thru 12 (representing 0.5 to 5.5 inches from the 
penetration point) are presented in Table IV and plotted ial 
Figs. 57 and 58. Table IV also shows the O, on the inside 


surface of the tank exit wall (wet) at the same time that 
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Maximum stress oF predicted by SATANS on 
the dry side of the plate plotted against 
distance from exit point. 
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the other stresses were measured. By comparing O, dry to 

©). wet, a profile of the stress through the plate thickness 
may be obtained. This profile reveals that significant 
stretching of the plate is predicted. Further examination 

of Table IV shows that the maximum stress predicted starts 

at station 2 and proceeds sequentially with time outward to 
station 12. Since the maximum stress around a crack is nor- 
mally at the crack tip, the predicted succession of maximum 
stresses moving outward from the center seems consistent with 
the idea of a crack growing from the penetration point outward 


toward the edge of the plate. 


B. PROCEDURE FOR PREDICTING THE MAXIMUM CRACK LENGTH 


Because the fracture surfaces of the actual test plates 
appear similar to simple tensile tests of thin plate specimens, 
tensile test data from the stretching of cracked plates was 
examined. For example, analysis of the data presented in 
Table CAD1O of Ref. 34 reveals a pattern in the fracture of 
0.125 inch thick 2024-T3 aluminum cracked plates. For test 
specimens of width 4, 8, and 24 inches, the stress required 
to cause a crack to increase in length was found to be 
a function only of the ratio of specimen width tO ergen 
length. The correspondence between the fracture stress and 
the ratio is shown in Table V and plotted in Fig. 59. 

Table VI shows the approximation made to the crack length 
face of Table Il in order to obtain a Single StraishtveEte aom 


equal damage to the small multiple cracks of the 7-HR series 
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TABLE V 


stress Required to Increase Crack 
Length in 2024-T3 Tensile Spee ae [34/7 


Specimen Width to 
Crack Length Ratio 


10/1 


10/3 
10/6 
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TABLE VI 


Length of Straight Crack Approximation of the 
Cracks from the 7-HR Series Low Velocity Shots 


APPROX. STRAICHE 


SHOT VELOCITY (fps) CRACK LENGTH (inches) 
7-HR-10 1347 4.0 
7-HR-11 1386 6.0 
7-HR-12 1366 Bar 
7-HR-14 1384 De 


low velocity snots. Since each test plate was 20 inches 

wide, the data of Table V gives a specimen width to crack 
length ratio of 10/2.5 for 7-HR-14. Referring to Fig. 59, a 
ratio of 10/2.5 corresponds to a stress of 35.6 ksi that is 
required to ae the 5-inch crack to increase. The 5-inch 
crack corresponds to a half crack length of 2.5 inches and 
station 6 of the computer model shown in Fig. 50. Since 

Ox acts perpendicular to the crack, the Ce stress at sta- 
tion 6 should be compared to the 35.6 ksi obtained above. 

The ©), Sireas of 11,3 kei predicted at Station 6 is jaeaeeron 
foes 5S Smaller than the 35.6 ksi. This correction factor of 
3.15 is, however, well within the range of correction factors 
obtained above by SOmoneTne predicted strains to actual strains, 
and therefore seems quite reasonable as a correction factor to 


be used in the prediction of the crack length. 


om 





Thus, the method of crack length prediction is as follows: 
the computer codes are used to predict the stresses in the 
exit wall. A correction factor, obtained ag explained above, 
is applied to the predicted stresses and the corrected stresses 
are plotted versus distance from projectile penetration eleva, « 
Assume various crack lengths and refer to Fig. 59 to find 
the stress required to increase crack length. The stresses 
required to increase the crack length are plotted on the 
same graph as the corrected predicted stresses. The pre- 
dicted crack length is the crack length at which the predicted 
stress first falls below the stress required for crack exten- 
Sion. For example, when the correction factor of 3.15 is 
applied to the O,.(dry) predictions listed in Table IV, the 
corrected stresses are those plotted in Fig. 60. Also plotted 
in Fig. 60 are the stresses required to increase the crack 
obtained from Fig. 59 for the 20-inch plate. Comparing the 
two stress plots reveals that the crack length will grow to 
5 inches, but no more than that because the stress required to 
extend the crack beyond 5 inches is greater than the predicted 
stress. The fact that the predicted crack length is equal to 
the actual crack length is of course, a consequence of the 
See Correction factor. 

Shots 10, 11, 12,and 14 were all similar enough to use the 
same computer stress predictions for all four shots. The 
proposed method of crack length prediction estimated the 


crack length to be approximately 5 inches long for each shot 
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as shown in Fig. 60. Examination of Table VI Shows that 
7-HR-11 would have exceeded the prediction by 1 inch. Howe 
ever, the other three shots had cracks that were less than 


the predicted crack length. 


C. EXAMINATION OF HIGH SPEED MOTION PICTURES 


High speed motion pictures of the exit wall response 
Were available for study for each of the NWC 7-HR series 
tests. The film speed was 8000 frames per second. From a 
close examination of the films it was possible to get an 
estimate of the lapsed time from initial penetration to the 
time when the cracks reached their maximum length. This 
data is presented in Table VII. All of the tests are not 
listed in Table VII because in some cases water spewing 
from the ruptured tank wall made it impossible to see the 
cracks. From Table VII it can be seen that the time it 
took for the cracks to grow to maximum length was between 
O.5 and 2.0 msec. The longer cracks taking the longer time 
to grow. As listed in Table VII, 7-HR-14 had a fracture 
time of 0.625 msec. Examination of Table VI shows that 
7-HR-14 has a crack approximately 5 inches long, or a half- 
crack length of 2.5 inches, which corresponds to station 6 
moebies, 5), SATANS predicted an alapeed time of 0.95 msec. 
for the maximum stress to reach station 6 after the projectile 
penetrated the plate. The 0.625 msec. observed on the films 


for 7-HR-14 compares quite well with the 0.95 msec. predicted 


by SATANS. 
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TABLE VII 


Estimated Lapse Time of Crack Process 
From High Speed Motion Pictures 


Estimate Lapse Time 


Snot Velocity of Crack Process (msec) 
7-HR-2 B27 1o25 
7-HR-3 2254 26 
7=HR-6 2763 1 ie 
7-HR-7 OS 256 
7-HR-8 2643 250 
7-HR-10 1347 Cr 
7-HR-12 1366 0.5 
7-HR-14 1384 0.625 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. FRACTURE SURFACE EXAMINATIONS 


The fractographic examinations of the damaged NWC 7-HR 
series test plates provided some information about the fracture 
process that took place as the cracks grew. However, it was 
not possible to come to any conclusions concerning the process 
or processes that took place at the point of projectile pene- 
tration. The projectile penetration point needs further study 
$n order to determine whether the plate fractures prior to the 
projectile penetration. 

The macroscopic and microscopic examinations of the frac- 
ture surfaces resulted in the conclusion that the fractures in 
the 7-HR series tests were generally ductile shear fractures, 
and were similar in many respects to the fractures that thin 
sheets in static tension tests exhibit. This similarity to 
tensile test specimens may be noteworthy since the time for 
the fracture process is generally under 2 msec. This may 
indicate that the strain rate effects of the fracture process 
are not important. Further investigation into the Straium 
rate effects are necessary to validate the use of static ten- 
sile test data in the fracture criterion. 

The fracture process itself deserves further investigation 
with different test parameters to confirm the characterization 


of the fractures as being ductile shear fractures and to more 


86 





firmly establish the similarity to static tensile test frac- 
tures. Eventually, temperature effects will also need to be 


considered because of the diverse temperature environment in 


which aircraft must operate. 


B. CRACK LENGTH PREDICTION METHOD 

This study proposes a macroscopic method that would allow 
the prediction of cracks caused by hydraulic ram phenonoma 
acting on the exit wall of a fluid contaanine tank. The 
method compares computer predicted stresses, which have been 
modified by an empirical correction factor, with tensile test 
@weeon the stress required to cause a crack in a pillate to 
increase. A detailed statement of the method is contained 
memoecetion V.B, The discussion in Section V.B. points out 
that the method provides accurate crack length predictions 
for the NWC 7-HR series low velocity shots. This close com- 
parison to the actual test results is very encouraging in 
view of the many simplifying assumptions that have been made 
in developing the method. Although the method predicts a 
crack of a certain number of inches, it should be used at 
the present time to obtain an idea of the relative size of 
the expected crack length as small, medium, or large: small 
being a crack length of about 25% of the plate width; medium 
being a crack length of between 25% to 50% of the plate width; 
and large being crack lengths larger than 50% of the plate 


width. 
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There are several areas that need further investigation 
in order to improve the method. First and probably most 
important the computer codes must be improved so that the 
predictions agree with the actual measurements and the cor- 
rection factor can be discarded. Also, the method needs to 
be applied to the medium and high velocity shots of the NWC 
7-HR series tests. Additional tests at each velocity level 
with different plate thicknesses are also required to validate 
the method. After the method has been validated for the 
2024-T3 aluminum considered in this study, other materials 


such as 7075-T6 aluminum must be examined. 


C. SIMPLIFIED DAMAGE PREDICTION 


In this study of the cracking of the exit wall a particu- 
lar feature has been noted. The damage to the exit walls 
caused by projectiles with initial velocity greater than 
approximately 2200 fps was considerably greater than the 
damage caused by projectiles at less than 2200 fps. This 
feature occurred on both the 0.125 inch thick 7-HR plates 
(Table III) and on the 0.063 inch thick 6-HR plates (Figs. 
2-5) tested by NWC. This appears to be more than coincidence 
and suggests that there is a consistent "ballistic velocity" 
that divides the damage into small and large categories. This 
feature should be investigated further since, if true, it 


would considerably simplify the analysis required to predict 


damage. 
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